Quenched QCD at finite temperature with chiral Fermions 
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We study physics at temperatures just above the QCD phase transition (Tc) using chiral (overlap) 
Fermions in the quenched approximation of lattice QCD. Exact zero modes of the overlap Dirac 
operator are localized and their frequency of occurrence drops with temperature. This is closely 
related to axial C/(l) symmetry, which remains broken up to 2Tc. After subtracting the effects of 
these zero modes, chiral symmetry is restored, as indicated by the behavior of the chiral condensate 
{{^ptp)). The pseudoscalar and vector screening masses are close to ideal gas values. 
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With new results from the Brookhaven heavy-ion col- 
lider appearing thick and fast Q , the time seems ripe for 
making a concerted effort to understand the dynamics of 
the high-temperature phase of quantum chromodynam- 
ics, namely the quark-gluon plasma. There are several 
puzzles that seem to have resisted a decade of efforts to 
understand them. The one we focus on involves the static 
screening of certain excitations of the plasma. 

It has long been understood that the screening of cur- 
rents in a plasma would give us information on its exci- 
tations. Currents with certain quantum numbers excite 
mesons from the vacuum at low temperatures, and should 
exhibit deconfinement related changes above the QCD 
phase transition temperature {Tc) ||5|]. Detailed studies 
have shown that this indeed does happen in the vec- 
tor, and axial- vector channels: the screening above Tc 
is clearly due to nearly non-interacting quark anti-quark 
pairs in the medium [|[0]. On the other hand, the scalar 
and pseudo-scalar screening masses show more compli- 
cated behavior — strong deviations from the ideal Fermi 
gas, and a strong temperature dependence. This puz- 
zling behavior is generic — it has been seen in quenched 
PI and dynamical simulations with two 19] and four fla- 
vors [p|- [7|,pO[ of staggered quarks, as well as with Wilson 
quarks [[ll[. This is the puzzle that we address and solve 
in this letter. 

The new technique we bring to bear on this prob- 
lem is to use a version of lattice Fermions called overlap 
Fermions [O. R has the advantage of preserving chi- 
ral symmetry on the lattice for any number of massless 
flavors of quarks jlj]. This is in contrast to other for- 
mulations such as Wilson Fermions which break all chi- 
ral symmetries or staggered Fermions which break them 
partially. Since the number of pions and their nature is 
intimately related to the actually realized chiral symme- 
try on the lattice, we should expect any realization of 
chiral Fermions on the lattice to provide insight into the 



question we address. 

The overlap Dirac operator (D) can be defined [Q in 
terms of the Wilson-Dirac operator (1?^,) by the relation 

D = l-D^,{DlD^)-^/^. (1) 

The computation oi D^^ needs a nested series of two 
matrix inversions for its evaluation (each step in the nu- 
merical inversion of D involves the inversion of Dl^D^). 
This squaring of effort makes a study of QCD with dy- 
namical overlap quarks very expensive. As a first step in 
this direction, we chose to work with quenched overlap 
quarks: to study the pattern of chiral symmetry restora- 
tion and screening masses at high temperature. 

We generated quenched QCD configurations at tem- 
peratures of T/Tc = 1.25, 1.5 and 2 on 4 x 8^ and 4 x 12^ 
lattices (see Table |). The corresponding couplings are re- 
spectively (3 = 5.8, 5.8941 and 6.0625. The configurations 
were separated by 1000 sweeps of a Cabibo-Marinari up- 
date. A previous computation |l5[| with T = quenched 
overlap quarks at the nearby couplings of /3 — 5.85 and 
6.0 allows us to compare finite and zero temperature 
physics. 

For the matrix M — Dl^D^, and a given source vec- 
tor 5, we computed y — M^^i'^h by a conjugate gradient 
(CG) version of a proposed Lanczos method ||l6[. CG 
gives better control over errors than Lanczos or other 
methods based on approximations using Chebychev or 
Reeves polynomials. In CG, the vector x = M^^b is ob- 
tained by iterating the residual vector r^-i-i = r^ — aiMpi 
and the direction of change Pi+i = Pi+iPi + ri+i, where 
cti = [ri,ri)/{pi,Mpi) and /J^+i = (r^+i, ri+i)/(ri, r^). 
The iterations are stopped when {ri,ri) < e, for a pre- 
determined tolerance e. Let Ns be the number of CG- 
iterations at stopping. In the orthonormal basis made of 
(li = '''i/ {fiifiY 1 the matrix M can be approximated 
by Q^TQ, where Q is made from the set {qi\ and T is 
a tridiagonal symmetric matrix of dimension N^. The 
non-zero elements of T are T^^j = l/cti -I- Pi/ai-i and 



7i^j_l_i = —yfPi/ai-i. Denoting by A the diagonal ma- 
trix of the eigenvalues of T, and by U its diagonalizer 
made from the corresponding eigenvectors, the desired 
solution isy = M-^/% ~ Q'^W K-^/'^UQb, when the CG 
iterations are started with ri:, = po = b. 




FIG. 1. Projection on the xy plane of two typical eigenvec- 
tors on a 4 X 12^ lattice at T = 1.5 Tc- The zero eigenvector 
(top) is strongly localized while the non-zero eigenvector (bot- 
tom) is not. 



A massive overlap operator is defined by 



D{ma) 



{l-ma/2)D, 



(2) 



where m is the bare quark mass, a the lattice spacing, 
and D is defined in p). We used the usual quark prop- 
agator, G{ma) = [1 - D/2]D-'^{ma) ^. We computed 
G on 12 point sources (3 colors and 4 spins) for 10 quark 
masses from 7na=0.001 to 0.5 using a multimass inver- 
sion of D'^D. The (negative) Wilson mass term in D^, 
which is an irrelevant regulator, was set to 1.8. The toler- 
ance was e = 10~^ in the inner CG and lO""* in the outer 
CG. A rough computation of the eigenvalues of D^'D, jjp' , 
was made on each configuration with a Lanczos method 
in each chiral sector. Whenever /i^ ~ O(10~^) was ob- 
tained, the few lowest eigenvalues and eigenvectors were 
refined to a precision of about 10~^ by a Ritz functional 
minimization. 

For each configuration we verified that the Ginsparg- 
Wilson ||l^ relation is satisfied to an accuracy of 10"^, 
thus ensuring a very precise implementation of chiral 
symmetry for our simulations. Another test of the preci- 
sion of our measurements was provided by a check of the 
chiral Ward identity — 
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where {'tpip) is the chiral condensate and xps is the pseu- 
doscalar susceptibility pq |. In all our computations we 
found that the equality was satisfied to better than 1 part 
in 10^. 

For most configurations we found that the spectrum 
of D'' D starts well away from zero. However, for some 



configurations we found zero and near-zero modes with 
fp ^ 10^'' clearly seperated by a gap from the non-zero 
modes with 0.1 ^ fi^ [|9[. The non-zero modes clearly 
came in degenerate pairs of opposite chiralities. For 
T ~ l.bTc and 2Tc the zero modes were all less than 10~^ 
and of definite chirality. Their numbers decreased rapidly 
with either increasing T or decreasing lattice size (see 
Table | for details). We constructed a gauge- invariant 
measure of localization pO] for a normalised eigenvector 
<i>(i) where i stands for position, spin and color, as the 
following sum. 
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This varies from unity for an eigenmode localized at just 
one site to 1/V for an eigenmode spread uniformly on 
a lattice of volume V. On 4 x 12^ lattices we found 
(T ~ (3 — 8) X 10"'^ for the zero modes and a < 10^'' for 
the non-zero modes (see Figure 0). At T = 1.25rc we 
also found near-zero modes with fj,^ < lO"**. These were 
as localised as the zero modes, but came paired in parity 
like the non-zero modes. 

The zero modes of chiral Fermions are related to 
instanton-like configurations ||2l| . These are, in turn, 
known to break axial U{1) symmetry. For two flavors, 
the order parameter for the axial C/(l) symmetry l22] 
is the difference of the flavor singlet and triplet scalar 
susceptibilities — 
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where n+ and n_ are the number of eigenvalues of D 
with left and right handed chiralities, and the difference 
needs to be evaluated only for the zero modes. In the 
quenched theory w, if non-zero, is singular in the ttj — > 
limit since ((n+ — n^Y') is independent of m. This is 
related to well known problems with x° in the quenched 
theory p3| . Combining our data with 124] , we find some 
evidence that ((n+ — ri-Y) /V falls as a high power of 
T/Tc- Nevertheless in quenched QCD, C/(l) symmetry is 
not completely restored even at 2Tc. 

Our measurement of ('0V') comes from the diagonal 
part of G{ma) on the 12 sources used for each configura- 
tion. The analysis proceeds by writing G in terms of the 
eigenvectors $^ of D'^ D with eigenvalue /i^ and chirality 
a, 



G. 



E^ 
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where A^ = 1— /i^/4. The contribution of a zero mode can 
be easily read off from the equation above and is seen to 
be proportional to $$ and l/ma. Since the eigenvector 



corresponding to the zero mode is localized, its contri- 
bution to the condensate depends strongly on the spa- 
tial position of the source vector. The remaining modes 
are delocalized and closely spaced; so the overlap of the 
eigenvector on the source is averaged out. 



0.84 

0.82 
I 0.80 - 
l|o.78 - 

V 

"to 0.76 
0.74 

0.72 - 
0.70 - 



FIG. 2. a^{Tpip)/ma as a function of ma at T — 1.25Tc 
(circle), l.STc (diamond) and 2Tc (square) on 4 x 12"^ lattices. 
Zero and near-zero mode contributions are subtracted. Data 
for 1.5Tc and 2Tc are displaced in ma for visibility. 

Our precise determination of the eigenvectors and 
eigenvalues allows us to subtract out the zero mode con- 
tributions in (g), although it can sometimes be a couple 
of orders of magnitude larger than the remainder. The 
subtracted condensate is strikingly identical to that seen 
in the sample without zero modes, at all the couplings 
and lattice sizes studied. We found that (tpip) varies lin- 
early with ma and goes to zero as ma -^ 0. This is how 
chiral symmetry restoration manifests itself in quenched 
QCD. 

In the thermodynamic and continuum limits, it is not 
clear whether the near-zero modes are related to chiral 
symmetry breaking; on any finite V they are not, but 
in these limits they may accumulate at zero. If they do, 
then (ipip) would not go to zero with m even above Tc 
in the quenched theory. Clearly, one needs to examine 
these limits very carefully. At T = 1.25Tc the number 
of near-zero modes is insignificant, but it seems that for 
T < 1.25Tc their numbers will be larger |gj| and the 
771 —> limit will have to be taken after the y — > oo 
limit. Quite likely, such a limit may have to be taken at 
more than one lattice spacing. 

For lattice Fermions which satisfy the Ginsparg- Wilson 
relation jl^, the following identities hold in the chirally 
symmetric phase as 77ia —f 0, 

Cs{z) ^ -Cpsiz) and Cv{z) = Cav{z). (7) 

Here C is the screening correlation function in the spa- 
tial z direction of an operator summed in the other three 
directions. The subscripts PS refer to a pseudo-scalar 
operator, 5* to a scalar, y to a vector and AV to an 
axial- vector [E5[. We find that the V and AV correlators 
indeed agree at all z and at all temperatures we stud- 
ied. This is in agreement with our conclusion that the 
temperature range we studied has chiral symmetry. 
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FIG. 3. The correlators Cps for configurations without 
zero modes (circle), Cps for configurations with zero modes 
(diamond) and {Cps — Cs)/2 for all configurations (square) 
for 4 X 12^ lattices with ma = 0.001 and T = l.STc. Some of 
the data points are displaced in z for visibility. 

It is clear from the chiral Ward identity (y) that large 
fluctuations in (ipip) due to the zero mode must also lead 
to similiar non-statistical fluctuations in Cps- Further, 
since a zero mode contributes identically to Cg and Cps , 
the two can even have the same sign if this contribution 
is large enough p6| . Ignoring the configurations with 
zero modes, the S and PS correlators do obey the iden- 
tity in (0), suggesting that simple results could be seen 
by eliminating the zero mode contribution. 

One way to do this is to use the measured $'s to 
subtract the zero modes and construct a new correlator 
Cps{z). Alternatively, we could consider the difference, 
{Cs{z) — Cps{z))/2, which should equal Cps{z). Figure 
H exhibits the comparison of the three correlators. Sim- 
ilar excellent agreement is also seen for these correlators 
for all 777a < 0.1 and at other T. 

After subtracting the effects of the zero mode, we find 
that the correlator identities are satisfied for the S/PS 
sector as well as V/AV. In addition, Cy is described well 
by an ideal gas of overlap quarks on the same lattice, as il- 
lustrated in Figure |^ (although the figure shows this only 
for one quark mass and T, this is true for all T > 1.25 Tc 
and all ma < 0.1). While Cps seems to differ from the 
ideal gas result, the measured values of the PS screening 
masses, Mps, are only 10% smaller than ideal gas screen- 
ing masses. This difference is small enough to be plau- 
sibly explained in a weak-coupling computation, quite 
unlike earlier results from staggered or Wilson quarks. 

Differences between T = mesons and our measure- 
ments are extremely clear. For T > 1.25Tc, Mps is con- 
stant and non- vanishing for 77ia ^0.1. In the same tem- 
perature range, the ratio mps{T)/mv{T) is within 10% 
of unity and quite different from the measured values at 
T = at nearby couplings /3 [|5|. Thus, the simple pic- 
ture that emerges is a property of the high temperature 
phase of (quenched) QCD. 




FIG. 4. The screening correlators on 4.12^ lattices at 
T — 1.5Tc for ma = 0.001, in analyses without zero modes. 
The V/AV correlators (boxes) agree very well with the ideal 
gas computation (solid line), and the S/PS correlators (cir- 
cles) are also similar. 

In conclusion, working with chiral (overlap) Fermions, 
we have found several new results and a consistent picture 
of the high temperature phase of quenched QCD. Axial 
C/(l) symmetry is not restored even at 2Tc. As a result 
the thermal ensemble contains gauge fields which give rise 
to Fermion zero modes of definite chirality. When the ef- 
fect of these modes is subtracted, {iptp) vanishes in the 
zero quark mass limit, showing that chiral symmetry is 
restored. Simultaneously, parity doubling is seen in the 
spectrum of screening masses, which are close to those 
expected in an ideal Fermi gas, even for the S/PS sector. 
Since some of these results are not obtained with stag- 
gered quarks, it is an interesting question whether the 
two flavor QCD phase transition is properly described 
by such a representation of quarks. 

Some interesting problems remain to be solved. At 
T < 1.25 Tc, there are near-zero modes. It cannot be 
ruled out that these modes shift the quenched chiral sym- 
metry restoration point away from Tc- However, this 
question is crucially related to the evolution of near-zero 
modes with lattice volume and spacing. Hence the nature 
of these complications, and the question of whether they 
are quenched artifacts or remain in full QCD, will only 
become clear with further studies which are underway. 
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Mps/T 
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1.25 
1.50 
2.00 


4 X 12'' 
4 X 12^ 
4 X 12^ 


50 
50 
50 


18 
8 
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4 




13.24 (9) 
13.31 (8) 
13.44 (6) 


6.0 (2) 
6.0 (1) 
6.0 (1) 


6.6 (1) 
6.5 (1) 
6.5 (1) 


1.50 
2.00 


4x8'' 
4x8^ 


100 
26 


1 







13.36 (5) 
13.34 (8) 


6.0 (2) 
6.0 (1) 


6.5 (2) 
6.5 (1) 



TABLE I. Temperatures (T), the lattice volume {V), and 
the number of configurations analysed (Nc). Also shown are 
the number of zero modes (No) and near-zero modes (Ag), 
and {ipip) /m, and the PS and V screening masses for m ^ 0. 



